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O Abstract 
(N 

Q We consider a quantum linear oscillator coupled at an arbitrary strength to a bath at an arbitrary 
(D 

Q temperature. We find an exact closed expression for the oscillator density operator. This state is 

O non-canonical but can be shown to be equivalent to that of an uncoupled linear oscillator at an 
effective temperature T*g with an effective mass and an effective spring constant. We derive an 

p I effective Clausius inequality dQ*g < T*g dS, where 6Q.*g is the heat exchanged between the effective 

^ (weakly coupled) oscillator and the bath, and S represents a thermal entropy of the effective oscil- 

d 

^ lator, being identical to the von-Neumann entropy of the coupled oscillator. Using this inequality 

(for a cyclic process in terms of a variation of the coupling strength) we confirm the validity of the 
second law. For a fixed coupling strength this inequality can also be tested for a process in terms 
fSl of a variation of either the oscillator mass or its spring constant. Then it is never violated. The 

in 

properly defined Clausius inequality is thus more robust than assumed previously. 

o 
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I. INTRODUCTION 



Thermodynamics of small quantum objects coupled to quantum environments in the low 
temperature regime has attracted considerable interest as the need for a better theoretical 
understanding increases in response to novel experimental manipulation of such systems. 
In particular, the finite couphng strength between system and environment gives rise to 
some quantum subtleties and so can no longer be neglected (calhng for methods addressed 
by 'quantum thermodynamics' [HISIIS]) whereas ordinary quantum statistical mechanics is 
intrinsically based on a vanishingly small coupling between them. 

At the heart of quantum thermodynamics, the foundational question as to the validity 
of the second law of thermodynamics comes up. In fact, with its challenge the applicability 
of thermodynamics is at stake. So far, the validity of this basic law has extensively been 
examined in the scheme of a quantum harmonic oscillator linearly coupled to an independent- 
oscillator model of a heat bath (quantum Brownian oscillator) in equilibrium at a low tem- 
perature T. It has been argued here that there is a violation of the Clausius inequality 
representing the second law [H [3] in such a way that SQ^TdS at T^O with respect to 
a variation of a Hamiltonian parameter of the coupled oscillator, namely, either its mass or 
spring constant. In the above relation, SQ is the heat and dS is the entropy change. 

However, following the second law in its Kelvin-Planck form, which states that it is 
impossible to devise a machine (i.e., a heat engine) which, operating in a cycle, produces no 
effect other than the extraction of heat from a thermal energy reservoir and the performance 
of an equal amount of work [6], it has been demonstrated that an apparent excess energy 
in the coupled oscillator at zero temperature (T = 0) is less than the minimum value of the 
work (equivalent to the Helmholtz free energy at a constant temperature) to couple the free 
oscillator to a bath so that the second law is not violated down to zero temperature [71 E]. 

This result has been generalized to a cychc process of couphng and decoupling between 
the oscillator and a bath at an arbitrary temperature by obtaining the positive-valuedness 
of the minimum work needed for the coupling minus the maximum useful work obtainable 
from the oscillator in the decoupling (the second law with respect to a variation of the cou- 
phng strength) |9]. This positive-valuedness is actually at its maximum at zero temperature 
and asymptotically vanishes with increasing temperature, whereas the classical counterpart 
would identically vanish at an arbitrary temperature (even for a non-vanishing coupling). It 
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was further claimed here that this quantum behavior is associated with the system-bath en- 
tanglement induced by the finite coupling strength between them (clearly, the coupled total 
system (i.e., the coupled oscillator plus bath) is in a thermal state with (partial) entangle- 
ment whereas the decoupled total system is simply in a separable state). It has, indeed, been 
found that at zero temperature the energy fiuctuation in the coupled oscillator can provide 
entanglement information [TOl E]. This claim was supported by the numerical analysis of 
the system-bath negativity as an exact entanglement measure |T2] that the negativity be- 
havior versus temperature is in accordance with the above quantum behavior of the second 
law up to the existence of the critical temperature above which the negativity vanishes. 

It has also been shown [I2] that the Clausius inequality (in terms of the equilibrium tem- 
perature of the total coupled system and the von-Neumann entropy of the coupled oscillator) 
is actually violated with respect to a variation of the mass of the coupled oscillator (not a 
variation of the coupling strength); the behavior of this violation versus the temperature is 
essentially different from that of the system-bath negativity so that it has been concluded 
that the system-bath entanglement is not responsible for the violation of the Clausius in- 
equality. However, as the reduced equilibrium density operator of the coupled oscillator is 
not in form of the canonical thermal state oc e~^^% there is not a well-defined /oca/ tem- 
perature of the coupled oscillator (especially in the low temperature hmit) so that applying 
the equihbrium temperature of the total coupled system for the violation of the Clausius 
inequality for the subsystem (actually with respect to a variation of a local parameter of the 
coupled oscillator, namely, either its mass or spring constant) is not justified. Further, this 
violation was actually based on the numerical findings [I2] that the heat 5Q exchanged with 
a bath in a reversible variation of the local parameter is always strictly greater than T dS, 
which, however, does not satisfy the equality condition of a well-defined Clausius inequality 
for the reversible process. 

On the other hand, introducing some generalized entropic measure and using its maximum 
condition [I3] it has been shown that the Clausius inequality obtained in some operational 
form is valid under such a generalization jH]. However, this approach is not directly appli- 
cable for the quantum Brownian oscillator since the reduced density operator of the coupled 



oscillator in equilibrium (c/. equations (17) and (18)) is not in form of the stationary state 
obtained from the maximum condition of this generalization. 

In this paper we intend to resolve the above controversial issue by introducing an effective 
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Clausius inequality with no violation, well-defined in the scheme of quantum Brownian 
oscillator. To do so, we begin with considering the reduced density operator of the coupled 
oscillator. 



II. REDUCED DENSITY OPERATOR OF THE COUPLED OSCILLATOR 

The quantum Brownian motion in consideration is described by the model Hamiltonian 

H = Hs + Hb + Hsb, (1) 



where 



"J 



N N 2 

j=i j=i ^ 

and the spring constants are ko = MuJq and kj = rrijUj'^. From the hermiticity of Hamilto- 
nian, the coupling constants Cj are obviously real-valued. The total system is assumed to 
be in the canonical thermal equilibrium state = e~^^ jZp where (3 = l/{kBT), and Zjj 
is the partition function. From the fluctuation-dissipation theorem [I5l[l6], it is known [T7] 
that 

= I j-duj coth(^) Im{x(^ + zO+)} (3) 
= ^ rduu' coth(^) Im{x(^ + zO+)} (4) 

in terms of the susceptibility 

TV 



n 

N 

n - ^i) 



xH = -^V (5) 



k=0 



where {uk} are the normal-mode frequencies of the total system H. For an uncoupled 
oscillator, Imx(t^ + «0+) obviously reduces to ^^o) and thus (g^)/? = coth 

and (p^)/3 = ^^^^ coth For the well-known Drude model (with a cut-off frequency ua 
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and a damping parameter 70), which is a prototype for physically realistic damping, we have 

m 

1=1 

1=1 

in terms of the digamma function ip{y) = d \nT{y)/dy [IB], where ui = 002 = zi, = Z2, 
and the coefficients 

_ ^1 + 22 . \(2) _ n + z2 . \(3) _ n + zi (o\ 
{n-zr){z2~n) ' "^d _o)(22-2i) ' ^d {z2 - n){zi - Z2) ' 

Here we have adopted, in place of (wq, u;^, 7o), the parameters (wq, fi, 7) through the relations 

m 



UJn := w 



■ u;,:=Q + ^ ■ 7o:=7 "?^7^r^ (9) 



and then zi = 7/2 + iwi and 2:2 = 7/2 — iwi with wi = a/ (wq)^ — (7/2)^. 

The equilibrium density operator of the coupled oscillator is known as [IH [I9] 



For an uncoupled oscillator this easily reduces to a well-known expression 

{q\pp\q') = tanh ^ exp (-f { (g + g')' tanh (^) + (g - q')' coth (^) }) (11) 
where the parameter 



Let us now derive a closed form of the matrix elements 

Pnm ■■= {n\ps\m) = j^^ j'^^dqdq' i)l{q) {q\ps\q') i)^{q') (13) 

in the basis composed of the eigenstates {\n),\m)} of an uncoupled oscillator to confirm its 
deviation from a (diagonal) form of the canonical thermal state pp. After making lengthy 
calculations, every single step of which is provided in a detail in Appendix we arrive at 
the closed expressions 

P2k,2l+1 = P2k+l,2l = (14a) 



in terms of the hypergeometric function 2Fi{a,b; c; z) [18]. Here, the four dimensionless 
quantities 



A, 



4c4 



(f^)^_ 1 



As shown, the reduced density matrix {ps)nm symmetric with respect to {n,m). Subse- 
quently we also find (c/. Appendix lAl) that for > / 



P2k,2l 



P2k+l,2l+l 



(-T^)' (A/i) (l-A,a) / r(fc+|) fc! r(2;+l) p(k-Lk-l) ( 1 

c^2m,A, Y T{i+\) n r(fc+m) ^21 l^^/I^ 



fc! r(2/+2) p(k-i,k-i) ( 1 

-to 



(16a) 
(16b) 



c^2(q^),A, Y r(i+|) r(fc+i+2) 21+1 

in terms of the Jacobi polynomial |IH] and for < / the matrix elements p2k,2i and P2fc+i,2Z+i 



correspond, respectively, to (16a) and (16b) with exchange of k and /. These are easily 



united into such a single expression that for either (n even > m even) or (n odd > m odd). 



r([=^] + i) [f ]! r(=±=.+i) ^™ 



n — m n — m ' 



(17) 

and for either (n even < m even) or (n odd < m odd) the matrix elements {ps)nm 



obviously given by (17) with exchange of m and n, where [y] represents the greatest integer 



less than or equal to y. On the other hand, from (14a) 



iPs)nn 







(18) 



where either (n even, m odd) or {n odd, m even). As seen, the reduced density matrix {ps)nm 
is, in general, not in diagonal form of the canonical thermal state e"'^^"*^""'"^)/^^ being valid 
for an uncoupled oscillator, where Ap q-P^o and A/j 0. This confirms that there 
is not a well-defined local temperature when one keeps starring at the oscillator Hg^ hence 



ignoring that it (strongly) couples to a bath (c/. Section III, in which, on the other hand, a 
well-defined effective local temperature is introduced). We note here, however, that by using 

/ra 1^ — 1) 



h 



2mu)o 



(a+a'), p 



-a), and the matrix elements (18) with a\n) 



and a)\n) = + 1 |n + 1), as is the case for an uncoupled oscillator, 

(v^Pn,n-l + Vn+ lpn,n+l) = 



Trs(gps 



h 

2mu)Q 



(19) 



n=0 
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and likewise (p)/3 = 0. Actually, it can straightforwardly be verified that = (p% = 
with I odd. 

Let us consider the probability of finding the nth eigenstate from the coupled oscillator 
in ps, which reads 



Pn — {Ps)nn 

in terms of the Legendre polynomial [IE 



P 4 



(20) 



n-2k 



(21) 



fc=0 



Here the normalization ^„Pn = 1 easily appears with the aid of (15) and the relation [T8 

oo 



n=0 



(22) 

where z = ^ and h = Apy^l — Ap. Then the internal energy of the coupled oscillator 



IS 



Us = {Hs), = Y.PnEn = %^ + ^ . (23) 

n 

A comment deserves here. In [lO] and then [21], each of diagonal elements Wn was obtained 



in terms of the Legendre polynomial P„ like in (20). On the other hand, the closed expression 



for off-diagonal elements {ps)nm (-'-'''' ^'^^ (-'-8) has not been known, while by means of 



the numerical integration of (13) the off-diagonal elements | {ps)nm I have been obtained for 



m, n < 10 in |2^. We will next consider the eigenvalues and eigenstates of the reduced 
density operator pg. 



III. EIGENVALUE PROBLEM FOR THE OSCILLATOR DENSITY OPERATOR 



The eigenvalue problem to be solved reads 

EV„ := /!lc?g'(g|ps|g')0n(g') = Pn4>n{q) 



(24) 



where the matrix element {q\psW) is given in (10), and the eigenvalue pn is the probability 
of finding the nth eigenstate of the coupled oscillator. Following the idea used in [22], we 
put an ansatz 

Mq') = .O^e-#'?^i7„(cg') (25) 
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[cf. (Al)) into the integral in (24), which will yield 

EV„ = 



exp 



2 + \n?{e)0 



g2 x(I). 



(26) 



Here f = y | + c^{q^)f3 + with v = \\/lJFhWTp^ and 

(I) := S^^dq'e-^^'-y)'H^{sci) 
where the dimensionless quantities 



q' 



q ; y 



The integral in (27) can be evaluated in closed form of [23 



(I) 



sy_ 



(27) 



(28) 



(29) 



For 0n(g) to be an eigenstate of the operator p^, from equations (24)-(26) with (29) we need 



to require the argument of the Hermite polynomial, -^j^p = —j){vy'v'^ — 2c^ {q'^)^) cq 
to equal cq, which immediately gives 

1/4 



and subsequently the probability for the nth eigenstate (pn{q) as 

Pn ■ 



1 \ n 



(30) 



(31) 



As a result, we obtained the eigenvalues pn and the eigenstates 4)n{q) in closed form. 



From comparison between (12) and (30), we introduce an effective mass Meg and an 



effective frequency Weg, which satisfy the relationship 



(32) 



For an uncoupled oscillator, this obviously reduces to {p^) p / {<f) p = {MujqY. For a later 
purpose, it is useful to note that either Meg or cjeff is not yet determined. The probability in 



(31 ) can then be rewritten as p„ = (1 — ■C/?)!^/?)" in terms of ^/j = {v — \)/{v + \) = e '^eff ^eff 



with an effective temperature 



and subsequently as 



A 



'eff 



P n -- z. 



1_ g-/3gff;!<^gff(n+|) 
eff 



(33) 
(34) 
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in terms of an effective partition function Z^s = e~^eff '«^eff("+5) = {csch(/?eff /ic<;eff/2)}/2. 
Therefore the density operator ps of the coupled oscillator is equivalent to that (peff) of an 
uncoupled linear oscillator 

at temperature T^s = l/{kB/3es) such that ps = e~^es^es /Zes. Here the spring constant 
kes = MeffCUgg. Needless to say, the temperature T^s — > T if the coupling constants cj 
in ([2]). With the aid of (32) and (33) we can easily confirm that 

which can subsequently be used for the internal energy 

[/eff = {Hes)p = ^ COth = CUeS^/WhWh (37) 

(note again, though, that since u^g is not yet determined, f/eg 7^ Ug). Further, choosing the 
effective frequency 



from (32), we then have U*^ = Us- Accordingly, M*^ = {p^)i3/Us and 

= l{ko + ^) . (39) 

Here /cq is the spring constant of the uncoupled oscillator Hg. All effective parameters are 
now uniquely determined in terms of the starred quantities, namely, M*^, k*^, the effective 
temperature T*g = —hujlg/{kB ln^/3), and the internal energy U*g = {H*g)[^{= Us) where 
H*ff = p^/(2Mgg) + (fcgg/2) (c/. note, on the other hand, that clearly ps = peS = pts)- With 
the aid of ^ and ([7]), figure [T] demonstrates that k^g > k^, which leads to M*^ > M from 



Mgg/Cgg = (j?) 13/ {(f) p. We will use the effective oscillator with (Mg^, k^^.T*^ in Section IV 
for a generalization of the Clausius inequality. 

Let us consider the thermal entropy of the effective uncoupled oscillator i^eff (of course, 
ffgg, too, as its special case), which is 

^eff {=Stg) = kslnZ^s + kBf3esUes = -A;^ {ln(l - ^^) + ^ In^/j} • (40) 

We can immediately verify that this is identical to the von-Neumann entropy of the coupled 
oscillator, 

Sn = kBiv + I) In [v+D- kB {v - I) In {v - I) (41) 
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[1] (note again that ps = PeS = Peff)- From figure [2| it is shown that Sn increases with the 
magnitude of the damping parameter and also with the temperature of the total system. 

Now we briefiy comment on the introduction of effective parameters: First, it has been 
shown in [22] that the coupled oscillator with (M, uq) at zero temperature of the total system 
(T = 0) can be interpreted as an uncoupled one with {M,ujes) in a thermal state with a 
finite effective temperature Teg (with (3 = oo), where tUeff = a/ {v^)oo/ {<f)oo/M. Using 
the very same technique, we have generalized this result into (Meg, cUeff; ^eff) in (|32)) and 



(33) for an arbitrary temperature of the total system. It is interesting to note here that 
f/eff ^ Usi(3 = oo) though, where U,s = {f)oo/i'2M) + {M(weff)V2}(g2)^. Secondly, it is 
also known [ElE^ that the coupled oscillator can exactly be seen as an uncoupled oscillator 
with an effective frequency 

u,s = ^ arccoth (I^I^MF)^) (42) 

and an eff'ective mass Meff = a/ {p'^) p / {q'^) p / uJeS in the canonical thermal state ps = e'^^" /Zp 
{i.e., Teff = T), which can be well understood simply as a special case of (Meg , tUeg, Teg) in 



(32) and (33) (note again that these are not the starred quantities). However, Uge 7^ Us 



either, whereas Us = U^g introduced below equation (38). 

It is also interesting to compare the internal energy U*g of the coupled oscillator with an 
alternative definition [71I81I91I251I261I27] 

W = -|lnZ^, (43) 

where the partition function = Tre~^^/Trbe~^^''. Here, Tr;, denotes the partial trace 
for the bath alone (in the absence of a couphng between system and bath, the function 
would exactly correspond to the partition function of the system only). Then it immediately 
follows that W = Us + {Hb)f3+{Hsb)p-{Hb)i3' ^ Us where (i^f,)^, = Tr^ [iib e-^^VTr, e"^^^) . 
Therefore the energy lA is not valid for the reduced system alone. The entropy S = 
kB(i.nZi3 — (3-^\iiZfs) can also be introduced here [261 121], which is, however, different 
from the von-Neumann entropy Sn (= S*^) for the reduced density matrix ps of the coupled 
oscillator. Actually, the entropy S cannot be derived from the Jaynes maximum entropy 
principle [28l [29] apphed for the reduced system whereas the entropy Sl^ can be so with 
the effective temperature T*^. As a result, all thermodynamic quantities resulting from the 
partition function Zj^ are not appropriate for the well-defined local thermodynamics of the 
reduced system. 
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IV. CLAUSIUS INEQUALITIES 



We will discuss the second law of thermodynamics in terms of the Clausius inequality. 



To do so, we need the relationship obtained from (23) 



dUs = J2 + ^^") (44) 

n 

where ^ Endpn = Trs{Hs dps) = SQs corresponds to an amount of heat added to the coupled 
oscillator, and YliPndEn = Trs{psdHs) = SWs an amount of work on the oscillator |^. For 
a later purpose, we first consider the well-defined Clausius inequality for a weakly coupled 
oscillator 

6Qs<TdS. (45) 

For a typical reversible process we have either a variation of the mass of the oscillator or a 
variation of its spring constant in such a way that 

MiQs = ^i,{p'), + f^e), (46) 

= ^(csch/^)^ = T^^ (47) 

ikQs = ^i-,{p'), + 'fi-^f), (48) 

= -ff (csch^)' = Tj^^, (49) 
respectively (c/. it is also noted that dWs/dM = -{p'^)p/{2M^) and dWs/dko = (g^)/3/2). 



We thus confirm the equality sign in equation (45). 

Next we consider the Clausius inequality in the process of couphng between oscillator 
and bath. The couphng process can be represented in terms of a variation of the damping 
parameter such that for a reversible process, 

f^Qs = f,U, = + f^if), (50) 



AS = fef^l„^ (51) 



where S = Sn and 
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Equation (50) can be evaluated in closed form with the aid of equations ([6])-([9]) such that 



f^,'' = it^^r (53) 



1=1 

3 

d'y 



-,{f)f = -MJ: {X'W + 2!£. A<« {-4 + J * {^)}'^) (54) 
where dui/d'y = and 



1=1 



dwg _ 1 i ^ 1 I _ 




(55) 



and 

-fir = {j^ + i * (^)} ^ + {-^ + s *'"(^)} If (66) 



in terms of the digamma function ^(y) and the trigamma function ip^^\y) = (P\nT(y)/dy 



2 



[T8] . Here we obtained equations (53) and ( [54j ) for the overdamped case (7/2 > wq), which 
is, still, found to hold for the underdamped case (7/2 < wq) as well, being expressed in terms 
of the functions with complex- valued arguments. Then we get a violation of the Clausius 
inequality, dQs/d'y > T 88^/9'^ as seen from figure |3] 

We, however, argue that this violation results from an inappropriate choice of temperature 
T being defined for the total system. We now propose a well-defined form of the Clausius 
inequality pertaining to the coupling process in such a way that 

5Ql.^ < T^^dS^ (57) 

where 50,*^ is the heat exchanged between the effective (weakly coupled) oscillator with 



(M*g, fcgg) in (39) and a bath at the equilibrium temperature T*^ = —hujl^g/ {kB\^ii3) with 
(38). For an reversible process we then have 



which can be shown to be identical to T^gdSN/d'y with the aid of (51) and (52). Therefore, 
there is no violation of the Clausius inequality! From U*g = Ug., we note that dU*^ = 
dUs = Us{'^) — Uq = AUs, in which Uq = {huJo/2) coth(/3/ic<Jo/2) for an uncoupled oscillator, 
and f dUs = § dU*g = 0, where <f = JJ + represents a cyclic process of the couphng and 
decoupling. From the first law dUs = SQs + dWs = 5Qlg + 6W*g with dWs/d'j = and thus 
d{Qlff - Qs)/d^ = -(9W*ff/(97, it also follows that 

f,Qs - i^ms = T:sIS^ (59) 
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where, with the aid of ([39|) 



eff 



Figure [4] shows that 5>Vgg/(97 < 0, which immediately leads to no violation of the inequal- 
ity dQs/d'y < T^gdSN/d'y. Here it should be noted that we have appropriately selected 



the effective oscillator H*^ with {M*g,u*g) from (32) to introduce an effective temperature 



Tgg without any ambiguity, which is now a critical element for the well-defined Clausius 



inequality in (57). 



Therefore, we are now in a position to understand, by means of the Clausius inequality 



(57), the validity of the second law in a cychc process of the couphng and decoupling between 
oscillator and bath at an equilibrium temperature T. The validity has actually been shown 
for zero temperature (T = 0) in [71 E] and later for an arbitrary temperature in j9] by 
verifying the second law in its Kelvin-Planck form ^j; it states that the minimum work AJF 
needed to couple the oscillator to a bath (in a reversible process), being equivalent to the 
Helmholtz free energy of the coupled total system minus the free energy of the uncoupled 
total system, cannot be less than the maximum useful work obtainable from the oscillator 
when it decouples from the bath such that 

AUs < AJ^ (61) 

(note the strict inequality and see below). Here we have on the left hand side the internal 
energy AUg as the maximum useful work obtainable from the oscillator on completion of 
the decoupling process. 

For the coupling-decoupling process (with a varying damping parameter 7' : — > 7 — > 0), 



inequality (57) can be transformed to ^ SQ*g/T*g < 0, which means, according to the Kelvin- 
Planck form, that the net work obtainable from the effective uncoupled oscillator (with an 
accordingly varying parameter 7') on completion of this cyclic process cannot be greater 
than zero. For a reversible process, this inequality then reduces to 

/;#j;^q:. + /:#?;^q:. = o, (62) 



which means that the minimum work {AUg) done onto the oscillator for Uq U*g on 
completion of the coupling exactly equals the maximum useful work releasable from the 

13 



oscillator on completion of the decoupling. For comparison, on the other hand, the free 
energy AJF is, by definition, the minimum work done on both oscillator and bath so that we 



get the inequality in (61 ). Anyhow, the second law holds in the coupling-decoupling process. 



Now we consider the Clausius inequality (with a fixed coupling strength) after completion 
of the coupling, which has been discussed so far, e.g., in [H O [12]. Note that we are now 
with the effective oscillator with {M*g, k*g) at temperature T*g. We can then show that for 
a reversible process. 



dM ^eff 
dko ^eS 



2 ^^ 



(63) 



2M*gdko\i^ IP ^ 2 dko\^ IP 



(41-^/3) (cschi^)' = T*,4^^,(64) 



which follow from (33), (36), (38) and (39), respectively. Therefore, there is no violation of 



the Clausius inequality at all! Note here as well that from the first law dUs = 6Qs + 5Ws 
5Qeff + "^^eff' effective work is also well-defined such as 



2M 



eff 



{dMMKn) {P')p + \ (dMMKs) {f)p 



(65) 



the closed form of which can immediately be obtained with the aid of M*q = {p'^)fj/Us and 



(39). 



For comparison, we next seek to have equations (46) and \A%\ in closed form in the Drude 



damping model. After making some lengthy calculations (c/. Appendix |B]), we then obtain 
the expressions 

3 



9MW 1 13 



1 lf,2\{d) , 1 T^{d) 



M 



1=1 



dM 



1 = 1 



dM J 



(66) 
(67) 



where 



and 



mJ^i ■- 

afco^y Ip 

dko \^ 1 13 



Puj,^ ^ 7T \ 2n J j dM + ^rf 



(0 



^ 2tt^ ^ y 2tt n dM 



du)i 



M X/ ^0^1 



(d) 



1=1 



-MY,{koK\''W + 2^A 
1=1 



2-K 



dko ) 



(68) 

(69) 
(70) 
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where koK^'^'' is defined as equation ([68j) but with the replacement ofd/dM by d/dko. Equa- 



tions (66) and (67) (as well as (69) and (70)) hold for both overdamped and underdamped 



cases. To discuss the second law, we now consider an equality similar to ( [59| in form of 



a 

aM 


Qs - 


a 

aM 


a 

ako 


Qs - 


a 

ako 



9 C 
^eff dM'^N 

rn-k a C 



(71) 
(72) 



obtained with the aid of the first law. Here dWs/dM = -{p'^)p/{2M^) and dWs/dko 



(g^)/3/2, whereas dWs/d'y for (59) vanishes. Actually we have here dM,ko(}Ves ~ > 0, 
which follows from dQs/dM > T^gdS^/dM and dQs/dko > T^gdS^/dko, demonstrated in 



figures [5] and [6| respectively (where equations (66)- (70) were used). This can be interpreted 
as follows. To define a well-defined (effective) local temperature of the oscillator, we need to 
"project" the coupled oscillator onto the effective oscillator. In doing so, it is required here to 
do additional work {W*g — Ws) onto the oscillator, whereas in the couphng process we need 
to release the work W*g from the oscillator. Without considering this work compensation 
we would consequently get a violation of the Clausius inequality. It is also interesting to 



rewrite (71) as 



a 

aM^s 



■Qs 



T^Sn + Y 



(73) 



in terms of the temperature T of the tota/ system, where Y = {T*g — T) dSN/dM + d{W*^ — 
Ws)/dM. Figure g shows that dQs/dM > TdSN/dM, which has been used in |12j for the 
justification of a violation of the Clausius inequality (note also the strict inequality even for 
a reversible process). However, we understand now that this simply represents the neglect of 
the additional term Y > rather than a violation proper. As a result, we have a generahzed 
form of the Clausius inequality, § 6Q*q/T*q < where § represents a cyclic process with 
respect to any variation of (M, ko, 7). For the relevant comment on effective thermodynamic 
relations, see Appendix [O 



V. CONCLUSION 



In conclusion, we have found a well-defined effective Clausius inequality appropriate for 
the quantum Brownian oscillator with any couphng strength. It satisfies the equality con- 
dition for a reversible process. We have clearly shown that there is no violation of the 
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inequality so that the second law of thermodynamics is robust even beyond the weak cou- 
phng limit. In doing so, we have used the effective internal energy U*g = {H*g)j3, being 
identical to the internal energy Ug = {Hs)i3, whereas the approach of apparently many other 



works has been based on a different energy U defined in (43) and discussed thereafter. 

We believe that this inequality will provide a useful starting point for a consistent gen- 
eralization of thermodynamics and information theory into the quantum and nano-system 
regime, respectively. As an example, a generalization of the Landauer principle |30l |3T] is 
in consideration, which can be understood as a simple logical consequence of the Clausius 
inequality; our findings then suggest an existence of an effective Landauer principle yet to be 
introduced rigorously |32] which is correct even in the strong couphng limit, whereas based 
on the violation of the Clausius inequality considered in [4, 5j as stated in Section |l| it was 
concluded in plj that the original form of Landauer principle may not be applicable in the 
strong coupling hmit. 

We thank A.E. Allahverdyan for his critical reading of the manuscript. One of us (LK.) 
is grateful to Prof. Jaewan Kim at KIAS for bringing Ref. [22] into his attention. We also 
appreciate all comments and constructive questions of the referees which made us clarify the 
paper and improve its quality. 



APPENDIX A: DERIVATION OF EQUATIONS (14a)-(14c) AND (16a)-(16b) 



Substituting into (13) the eigenfunction 



Mq) = ^2^6-^" H^icq) (Al) 
in terms of the Hermite polynomial where u = n,m, then we immediately obtain 



^ ^ JZIZodydy'H4y)HUy')x 



exp {-a,(y2 + y-) + ^ - ^) yy'} (A2) 

where y = cq and y' = cq', and ap = ^ + ^^^2^^ + The substitution of the relations 

m 

Hn{y) = (I)" e2^*-*%=o (A3a) 
HM) = ar e'y'^-^'\s=o (A3b) 
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into (|A2| subsequently allows us to have 



Pn 



(A4) 



where 



2c2 U(g2)^ 



?i2 



s. By using the identity 



J —00 



e 



(A5) 



we can first carry out the integration over y' in (A4) and then over y, which will give rise to 



1 1 



\/2"+'" n! m! c ^2(g2)^A^ 

/ (Y^)'" " (A^)" ^ .Q_.n 

V 2"+™n!m! c ./^T^^T^ 



t,s=0 



u=0 



in terms of A/s, Ta, Aa and in (15). Here we used equation (A3b), and u 



(A6) 



From the Heisenberg uncertainty relation with {q)f3 = {p)[3 = (cf. (19)), it follows that 
1 > A/j > 0. For a later purpose it is useful to confirm that for an uncoupled oscillator, 



= = and A 



13 = e 



To arrive at a closed form for pnm, we consider the expression in (A6) 



nm — 'P 



n 



r=0 



u=0 
(n-r) 



{Hm{u)} 



(r) 



(A7) 



u=0 



where (") = (^J^,^^, , and (■ ■ ■ = (|^)'' (■■■). The Hermite polynomial Hm{u) = Ji^u'^ + 
mhm-2 u"^"^ + ■ ■ ■ can be expressed as [T8] 

(A8a) 
(A8b) 

in terms of the confluent hypergeometric function iFi{a;b;z) = YltLo r(fe+fcj ^ with 
0, /c = 0, 1, 2, ■ ■ ■ . Then it immediately follows that 



r(-fc) 



«=0 



{2p) \ 2ih2p for r = 2p even, where p < 
otherwise 



(A9) 
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and from (|A8al) 



W{20! (2p)! (-0, 



where the Pochhammer symbol {z)k = ^^i^^- Similarly, we can also find that 



u=0 



{2p + 1)! 2i+ih2p+i for r = 2p + 1 odd, where p < I 
otherwise 



and from (ASb) 



u=0 



{21+1)1 i2p+iy. (-Op 

^ ^ n p! r(p+|) ■ 



Also, in (A7) we have 



(n— r) 



(-A;3)^(2g)! 



for n — r = 2q even 



u=0 



otherwise. 



With the aid of (A9)-(A13), equation (A7) reduces to 



n,2l^f3 



k ( —1 

v^(-l)fe+' (2fc)! (20! (A^)'-' x-^ r/^i^„\ (a 



= (-1) 

for n = 2k even, and similarly 



n r(-o 



k+l (2k)\ (20! 

k\ v. 



p=0 



(|+p)r(fc-p+i) p! 



n,2Z--/3 







-o/- iNfe+i (2fe+l)! (2Z-M0! 
n,2;+l--/3 - ^ l-ij — fc! 



for n = 2A; + 1 odd. In (A15) we used the identity 



(AlO) 



(All) 



(A12) 



(A13) 



(A14) 
(A15) 

(A16) 
(A17) 



±il^(A,)SFi(-fc,-/;|;^) (A18) 



1 _ (-i)pr(-fc+p) /A.qN 
r(fc-p+i) r(fc+i)r(-fc) l^-^^J 

r(c) Y^oo r(a+fc)r(fc+fc) ^fc 



to get the hypergeometric function 2Fi{a,b] c; z) = r(a)r(fe) XlfcLo r( 



c+fc) fc! 



From (A6), (AT) and (A14)-(A18) with the relation jl8 



T{2u) = ^22-ir(z.)r(z.+ i) 



(A20) 



where u = k,l, we finally get equations (14a)-(14c). We will below simplify the closed forms 



in (14b) and (14c), respectively. 
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Now we use the relation [18] 
to express the matrix elements {ps)nm terms of the Jacobi polynomial [T8 

n 



k=0 



where /x, z/ > —1. Equation ( |A21[ ) allows us to have 

-k 



2F1 ( -k, -I; \] ^ 
2FA-k,-l 



■ 3. _L 

' 2' 



r(|)r(f+i) 
r(|)r(i+i) 



fovk>l where v = 1 — j^^- Further, we find that [33] 



^0 



, r(2n+i.+l)r(n+l) p(-5.^)/-| 



-1 



Xn+i/+l)r(2n+l) " 
,„ r(2n+iy+2)r(n+l) ^ pd-'") 



r(n+i/+l) r(2n+2) 



zP^5"'^(l-2z^ 



(A21) 
(A22) 

(A23a) 
(A23b) 

(A24a) 
(A24b) 



which can be verified, respectively, by using equation ( |A22[ ) and then comparing each coeffi 



cient oiz^ on both sides. Note here that P'^^^^^^'yz) = P. 



2n 



-z) andPi,'^-'^)' 



With the aid of (A23a)-( A24b), equations (14b) and (14c) can then be transformed into (16a) 



and (16b), respectively. 



APPENDIX B: DERIVATION OF EQUATIONS (66)-(70) 



Using the relations in ([o]) with duo/dM = —ujq/{2M) we can easily obtain 



Wo7 



. dj_ _ _dn_ . dwo _ wor2(wg+0^-7^) 

dM ' dM 2M(f2+7){n(7-f2)-w2} 



dM n+7M{n(7-n)-w2} ' dM 

From dujo/dko = —{duoQ/dM)/uo1, it immediately follows as well that 



fBll 



do. 
dko 



1 n+7 m_ . dj_ 

(wo)2 a dM ' dko 



dko ' dko 



1 n+7 awo 

(wo)2 n 9M 



We also have dui/dM = dVt/dM and 



(B2) 



d u!2 
dM 



d u)-i 
dM 



dwQ 



' ^ \2 dM 



dn 
dM 



9wo 1 



s 2 9A/ 
-1 



1 + 



' 2wn \2 



(B3) 
(B4) 
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And from ([s]), 



dX 



(1) 



and dX^f^ /dM —>■ dX]^' /dM with (2:1 ^ ^2)- Then we finally arrive at the expressions in 



{Q--2i)2 (n-z2)2 



.(2) 



' dM 



(21-22)2 



(B5) 



(66) and mjn, respectively. With the replacement oid/dM in (B3)-(B5) by d/dko, we can 



also have (69) and (70). 



APPENDIX C: COMMENT ON EFFECTIVE THERMODYNAMIC RELATIONS 

In ordinary thermodynamics, the notion of temperature T appears conceptually as a 
partial derivative of internal energy U with respect to entropy 5* such that T = dU/dS. This 
also holds for the effective temperature T^g-. Combining the first law dUg = SQ^g + 6W*g 
with the second law SQ^g = T*gdSN for a reversible process, we have 



^eff 



dUs 



It is also interesting to note that equation (37) can be recovered by the relation 



(CI) 



Us = utg 



(C2) 



which is well-defined in terms of a generating function 



eff 



eff 



(C3) 



From (40) it then follows as well that the effective free energy F^g = Us — T^gSjy, where 



eff 



-keT^glnZ. 



eff- 
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Fig. |T] (Color online) y = ko/k*e = 2/(1 + (^^ + 7) {p^)f^/{^ {Mwof {q^)p}) versus 
X = kBT/hwo (dimensionless temperature), where wq is the renormahzed eigen frequency 



of the oscillator Hs (cf. equation (jOj)); for k*^ refer to equation (39). From bottom to top, 
(blue sohd: damping parameter 7 = 10, overdamped), (black dash: 7 = 4, overdamped), 
(green solid: 7 = 3/2, underdamped) and (red dash: 7 = 1/2, underdamped). We have 
y<lso that M/M*g = 1/(2 - y) < 1. Here h = ke = wq = n = M = 1. 



Fig. [2| (Color online) y = Sn (von-Neumann entropy) versus x = kBT/hwo (dimension- 
less temperature); for Sn refer to equation ( [4l| . From top to bottom, (blue solid: 7 = 10), 
(black dash: 7 = 4), (green sohd: 7 = 3/2) and (red dash: 7 = 1/2). As 7 decreases, then 
Sn decreases. Here /i = /c^ = wq = f2 = 1. 



Fig. |3j (Color onhne) y = 10 ■ {dQs/d'j — T OSn /d'j) /h (dimensionless) versus x 



kBT/hwQ (dimensionless temperature); for y refer to equations (50) and (51). From bottom 
to top (at T = 0), (blue solid: 7 = 10), (black dash: 7 = 4), (green sohd: 7 = 3/2) and (red 
dash: 7 = 1/2). Here h = ks = wq = = M = 1. 



Fig. 14} (Color onhne) y = 100 ■ (dW^g/d'y) /h (dimensionless) versus x = ksT/hwo 



(dimensionless temperature); for y refer to equation (59). From top to bottom (at T = 0.5), 
(blue dash: 7 = 10), (red sohd: 7 = 1/2), (black dash: 7 = 4) and (green solid: 7 = 3/2). 
Here h = ks = wq = Q = M = 1. 



Fig. [5} (Color online) y = {dQs/dM — T*gdSN/dM) /{hwo/M) (dimensionless) versus 
X = ksT/hwo (dimensionless temperature); for y refer to equation ( [7l| ). From top to bottom, 
(blue sohd: 7 = 10), (black dash: 7 = 4), (green solid: 7 = 3/2) and (red dash: 7 = 1/2). 
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As 7 decreases, then y decreases. Here h = ks = = ^ = M = 1. 



Fig. [6j (Color online) y = {dQs/dko — T*gdSN/dko) / {h/ {Mwq)) (dimensionless) versus 
X = ksT /hwQ (dimensionless temperature); for y refer to equation ( [72) ). From top to bottom, 
(blue solid: 7 = 10), (black dash: 7 = 4), (green solid: 7 = 3/2) and (red dash: 7 = 1/2). 
As 7 decreases, then y decreases. Here h = ks = = ^ = M = 1. 



Fig. [7| (Color online) y = [dQ^/dM — T dS^/dM) /{hwo/M) (dimensionless) versus 
X = ksT /Hwq (dimensionless temperature); for y refer to equation (73). From top to 
bottom, (blue solid: 7 = 10), (black dash: 7 = 4), (green solid: 7 = 3/2) and (red dash: 
7 = 1/2). As 7 decreases, then y decreases. Here h = kB = wq = = M = 1. 
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